Improving imaging resolution of shaking targets by Fourier-transform ghost diffraction 



O 
<N 



Cong Zhang, Wenlin Gongj3 and Shensheng Han 
Key Laboratory for Quantum Optics and Center for Cold Atom Physics of CAS, 
Shanghai Institute of Optics and Fine Mechanics, 
Chinese Academy of Sciences, Shanghai 201800, China 
(Dated: July 27, 2012) 

For conventional imaging, shaking of the imaging system or the target leads to the degradation 
of imaging resolution. In this work, the influence of the target's shaking to fourier-transform ghost 
diffraction (FGD) is investigated. The analytical results, which are backed up by numerical simula- 
tion and experiments, demonstrate that the quiver of target has no effect on the resolution of FGD, 
thus the target's imaging with high spatial resolution can be always achieved by phase-retrieval 
method from the FGD patterns. This approach can be applied in high-precision imaging systems, 
to overcome the influence of the system's shaking to imaging resolution. 
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Relative motion between the imaging system and the 
object is one of the common reasons for imaging degra- 
dation. For instance, camera's shaking or object's high- 
speed moving would cause the blur and degradation of 
imaging resolution. For conventional imaging, a usual 
method to decrease motion blur is to reduce the expo- 
sure time, while this will reduce the signal-to-noise ratio 
(SNR) of the images. 

Other approaches for motion deblurring include me- 
chanical compensation, optical compensation and elec- 
trical compensation Q. However, the imaging reso- 
lution is limited by the compensation accuracy. Also, 
the designing and manufacturing of compensation de- 
vices are complex and the development of compensation 
techniques is still one of the most challenging tasks. 

Motion deblurring can also be completed by digital 
image processing method [3]. The blur kernel, namely 
point spread function, is the key factor in the process of 
image restoration. If the blur kernel is known before im- 
age restoration, inverse filtering and Wiener filtering are 
widely used for image restoration [4-6] . If the blur kernel 
cannot be known or measured in some practical applica- 
tions, image restoration is usually realized by blind de- 
convolution technique (^-Q- However, image restoration 
needs large computation and suffers from serious restora- 
tion errors, moreover it is sensitive to noise disturbance. 

Ghost imaging (GI), as a nonlocal imaging method, 
has been widely investigated during last fifteen years [10- 
fl7| . An interesting and obvious feature of ghost imaging 
is that fourier-transform imaging of the object can be ob- 
tained by measuring the intensity correlation function be- 
tween two light fields, even when the source is a thermal 
light source and the object is located in Fresnel region 
relative to the source [l7l-[22j. Recently, super-resolution 
far-field ghost imaging and ghost imaging lidar are also 
obtained by exploiting the sparsity constraints of images 
(23l-[25|. However, all the previous work on ghost imag- 



ing is focused on the static targets. Obviously, imaging 
an unstable and moving target is meaningful in practice. 
Based on the lensless fourier-transform ghost diffraction 
(FGD) scheme [20] and phase-retrieval techniques [2lj, 
the effect of the target's shaking on FGD is investigated 
both theoretically and experimentally, and the target's 
imaging is retrieved from the experimental FGD pat- 
terns. We show that the scheme is helpful in overcoming 
the shaking blur due to the novel characteristic of ghost 
imaging technique. 
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FIG. 1: Schematic of lensless FGD for shaking targets with 
pseudo- thermal light. The object is operated by a motion 
platform so that it can shake on the object plane, perpendic- 
ular to the optical axis. 



Fig. 1 presents the experimental schematic of lensless 
FGD for shaking targets with pseudo-thermal light. The 
pseudo-thermal source, which is obtained by modulating 
a laser beam (the wavelength A=532 nm and the light 
spot's transverse size 1^=4 mm) with a rotating ground 
glass, is divided by a beam splitter (BS) into a test and 
a reference paths. In the reference path, the light propa- 
gates directly to a charge-coupled device (CCD) camera 
Dr. In the test path, the light goes through an object 
and then to a single pointlike detector Dt. The object is 
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placed on a motion platform and the platform is driven 
by a stepping motor so that it can move straight in one 
dimension. In addition, the stepping accuracy of the plat- 
form is 5 jiim and the stepping motor is connected to a 



computer, which can control the motion of the object. 

Based on GI theory [l5|, [I^, the correlation function 
between the two detectors is 



AG^^^^\xr,Xt) = {Er{Xr)E:{Xr)E:{Xt)Et{Xt)) - {Er{Xr)E:{Xr)) {Et{Xt)E: {Xt)) 

= J dxoidxo2dxosdxo4 (£^0(^01)^^(^04)) (£^0(^02)^0(^03)) 

X {hr{Xr,X0i)hl{xt,X04)K{^r,X02)ht{xt,X03)) , 

I 



(1) 



where < • • • > denotes an ensemble average, Er{xr) and 
Et{xt) denote the light fields on the reference and test de- 
tection planes, respectively. {Eq{xo)Eq{xq)) is the auto- 
correlation function of light field on the ground glass 
plane. ht{xt,xo) is the impulse response function of the 
test path, whereas h^{xr^ xq) denotes the phase conjugate 
of the impulse response function in the reference path. 

For the schematic shown in Fig. 1, under the parax- 
ial approximation, the impulse response function of the 
reference path is 



(2) 



where k = 27r/A. And the impulse response function for 
the test path is 



ht{xt,xo) (X [ dx't{x' ,T)ey.-^{^{x' - xqY) 
J 

X exp{|— (xt -X )^}, 



(3) 



where t{x\r) is the transmission function of the object 
at time r. Assuming that the light field on the ground 
glass plane is fully spatially incoherent and the intensity 
distribution is uniform as a constant intensity Iq^ then 



{Eo{xoi)E^{xo4)) = IoS{xoi - X04), 

(£^0(^02)^^0 (^03)) = IoS{xo2 - ^03), 



(4) 



where 6{x) is the Dirac delta function. 

Substituting Eq. (2)- (4) into Eq. (1) and suppose that 
z = zi -\- Z2^ the correlation function is 

AG*^^'^^(xr,Xt) (X / dx[dx2exp{ — — {xr — xt)x[} 

J ^2 

X (t(xi,r)t*(x2,r))exp{ — (x^ - Xt)x2}, (5) 

Z2 

The object's shaking means that the center position 
of the object will deviate from the optical axis. If the 
object's center position deviates ^ from the optical axis 
at time r, then the object's transmission function t{x' ^ r) 
will change into t{x' — C)- practical applications, al- 
though the shaking deviation of the object is unknown 



and random at every time, it obeys to a certain statistical 
distribution. Assume that the shaking deviation ^ obeys 
to probability distribution P(^) (where J (i^P(^) = 1), 
then 

{t{x[,T)t*{x'^,T)) = J dCP{0t{x[-0t*{x'^-0, (6) 

Substituting Eq. (6) into Eq. (5), after some calcula- 
tions, Eq. (5) can be represented as 



AG^2'2)(x„xt) oc 



AZ2 

f{^{Xr - Xt)) 



(7) 



where T{q) is Fourier transformation of the object t{x'). 
If the detector in the test path is a single pointlike de- 
tector which is located at = 0, then 



AG^^^'^\xr,Xt=0) (X 



(8) 



From Eqs. (7) and (8), the object's FGD can be still 
obtained even if the object is shaking and the result is the 
same as that achieved in the static condition [17]. There- 
fore, the object's shaking has no effect on the imaging 
resolution of FGD. 

In order to verify the analytical results. Fig. 2 and 
Fig. 3, respectively, present the numerical simulated and 
experimental results of imaging a three-slit in different 
shaking conditions, using the lensless scheme shown in 
Fig. 1. In the simulation and experiments, the distances 
listed in Fig. 1 are as follows: z = 550 mm, zi = 200 
mm, Z2 = 350 mm. The object is a three-slit with slit 
width a = 0.2 mm, center-to-center separation d = 1.0 
mm, and slit height h = 1.0 mm. 

The shaking modes of the object are controlled by the 
motion platform and two statistical distributions are dis- 
cussed: uniform and normal. The maximum shaking de- 
viation from the optical axis is 100/im, 500/im, 1000/im, 
2000/im, respectively, which is used to perform the infiu- 
ence of the shaking amplitude to the imaging resolution 
of FGD. 
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FIG. 2: Numerical simulation results of shaking object FGD 
with the lensless scheme shown in Fig. 1 (averaged 4000 mea- 
surements). X-axis indicates the position on the reference 
detector Dr. Y-axis indicates FGD reconstructed by mea- 
suring the intensity correlation function AG*^^'^-* (xr, xt = 0). 
(a)-(d): the shaking mode obeys to uniform statistical dis- 
tribution, and the maximum shaking deviation from the op- 
tical axis is lOO/xm, 500/im, lOOO/xm, 2000/im, respectively. 
(e)-(h): the shaking mode obeys to normal statistical distri- 
bution, and the maximum shaking deviation from the optical 
axis is respectively lOO/xm, 500/im, lOOO/xm, 2000/im. Red 
dashed curves present FGD reconstructed in different shaking 
conditions versus black solid curves for the static condition. 



For Fig. 2 and Fig. 3, FGD in the static condition is 
illustrated by black solid curves. If the object's shaking 
mode obeys to uniform distribution, the results in differ- 
ent shaking amplitudes are shown in Fig. 2(a)-(d) and 
Fig. 3(a)-(d). While in the case of normal distribution, 
Fig. 2(e)-(h) and Fig. 3(e)-(h) present the corresponding 
simulated and experimental results. From the simulated 
and experimental results shown in Fig. 2 and Fig. 3, 
both the shaking mode and the shaking amplitudes of 
the object have no effect on the imaging resolution of 
FGD and the imaging quality is the same as that in the 
static case, which is in accordance with the analytical 
results described by Eqs. (7) and (8). 



FIG. 3: Experimental results of shaking object FGD, the con- 
ditions are the same as Fig. 2. The images in different shaking 
conditions reconstructed by measuring the intensity correla- 
tion function AG^'^''^\xr, xt — 0) are listed at the upper left 
corners and red dashed curves are their cross-sections. In con- 
trast, black solid curves are the cross-section of FGD obtained 
in the case of static condition. 



In addition, based on phase-retrieval iterative algo- 
rithm [21], the object's imaging in two shaking modes 
and different shaking amplitudes, as shown in Fig. 4, 
can be retrieved from the experimental FGD patterns of 
Fig. 3. From Fig. 4(d,f) and (i,j), even if the object's 
shaking amplitude is much larger than center-to-center 
separation of the object, the spatial resolution of recon- 
structed object is nearly the same as that retrieved in the 
static condition when FGD technique and phase-retrieval 
method are combined, which can overcome the influence 
of relative motion between the imaging system and the 
object to imaging resolution in conventional imaging sys- 
tem. 

For FGD, a single pointlike detector far from the ob- 
ject in the test path can receive the global information 
of the object. However, the spectrum center of FGD will 
move if we change the position Xt of the single pointlike 
detector, which can be explained by Eq. (7) because the 
center of the object's FGD is related to the detection po- 
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FIG. 4: The object's imaging retrieved from FGD patterns 
based on the experimental results of Fig. 3. (a) The orig- 
inal object used in numerical simulation and experiments; 
(b) phase-retrieval image of static object FGD; (c)-(f): cor- 
responding phase-retrieval results of Fig. 4(a)- (d); (g)-(j): 
corresponding phase-retrieval results of Fig. 4(e)- (h). 



sition xt and the result has also been demonstrated in 
Ref. [26j. Different from the case of changing the detec- 
tion position Xt^ changing the object's center position has 
no effect on FGD. As a result, the technique of FGD can 
overcome the influence of relative motion between the 
imaging system and the object to the imaging resolution 
and the object's imaging with high spatial resolution can 
be always reconstructed by phase-retrieval method from 
the FGD patterns, which provides a brand-new approach 
of motion deblurring and is very useful to high-precision 
imaging systems and imaging of moving target. 



In conclusion, we have given a proposal to improve the 
imaging resolution of shaking object using ghost imag- 
ing technique. Both analytical and experimental results 
have shown that the object's shaking has no effect on the 
imaging resolution of FGD and the object's imaging with 
high spatial resolution can be always reconstructed from 
the FGD patterns. This technique has a great applica- 
tion in practice to overcome motion blur. 
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